Surface density of states and topological edge states in non-centrosymmetric 

superconductors 
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We study Andreev bound state (ABS) and surface density of state (SDOS) of non-centrosymmetric 
superconductor where spin-singlet d-wave pairing mixes with spin-triplet p (or /)-wave one by spin- 
orbit coupling. For d xy + p-wsve pairing, ABS appears as a zero energy state. The present ABS is 
a Majorana edge mode preserving the time reversal symmetry. We calculate topological invariant 
number and discuss the relevance to a single Majorana edge mode. In the presence of the Majorana 
edge mode, the SDOS depends strongly on the direction of the Zeeman field. 

PACS numbers: 74.45.+C, 74.50.+r, 74.20.Rp 



INTRODUCTION 

Recently, physics of non-centrosymmetric (NCS) su- 
perconductors is one of the important issues in condensed 
matter physics. Actually, several NCS superconductors 
have been discovered such as CePtsSi [if, Li 2 Pt3B Q 
and LaNiC2 0- Also, the two-dimensional NCS super- 
conductivity is expected at the interfaces and/or surfaces 
due to the strong potential gradient. An interesting ex- 
ample is the superconductivity at LaAlOs/SrTiOa het- 
erointerface Q. In NCS superconductors, the spin-orbit 
coupling comes into play. One of the remarkable features 
is that due to the broken inversion symmetry, supercon- 
ducting pair potential becomes a mixture of spin-singlet 
even- parity and spin-triplet odd-parity Frigeri et al. 
@ have shown that p (p x ±ip y )-w&ve pairing state has the 
highest T c within the triplet-channel in CePtsSi. Due to 
the mixture of singlet s-wave and triplet p-wave pairings, 
several novel properties such as the large upper critical 
field are expected @, Q • 

Up to now, there have been several studies about su- 
perconducting profiles for NCS superconductors 
In these work, pairing symmetry of NCS superconduc- 
tors has been mainly assumed to be s + p-wave where 
spin-triplet p(p x ± ipj,)-wave and spin-singlet s-wave pair 
potential mixes each other as a bulk state. However, in 
a strongly correlated system, different types of pairing 
symmetries are possible. Microscopic calculations have 
shown that spin-singlet d x 2_ y 2-w&ve pairing mixes with 
spin-triplet /-wave pairing based on the Hubbard model 



broken inversion symmetry. Also, a possible pairing sym- 
metry of superconductivity generated at heterointerface 
LaAlC>3/SrTi03 |4| has been studied based on a similar 
model [15| . In Ref. (ljj . it has been found that the gap 
function consists of spin-singlet d xy -w&ve component and 
spin-triplet p-wave one. The ratio of the d xy -w&ve and 
the p x (p y )-w&ve component in this d xy + p-wavc model 
continuously changes with the carrier concentration. 

Stimulated by these backgrounds, a study of Andreev 
bound states of d xy + p or d x 2_ y 2 + /-wave pairing has 
started 16j. It has been known that the generation of 



Andreev bound state (ABS) at the surface or interface is 
a remarkable feature specific to unconventional pairing 
since ABS directly manifests itself in the tunneling 
spectroscopy. Actually, for d xy -wave pairing, zero energy 
ABS appears [Hj]. The presence of the ABS has been ver- 
ified by tunneling experiments of high-T c cuprate IH 19[ 
as a zero bias conductance peak [20] . For chiral p-wave 
superconducting state realized in Sr 2 Ru04 [2l|, ABS is 



generated as a chiral edge mode which has a dispersion 
proportional to the momentum parallel to the interface 
[22J. For s + p- wave NCS superconductors, when the 
magnitude of p-wave pair potential is larger than that of 
s-wave one, it has been shown that ABS is generated at 
its edge as helical edge modes similar to those in quantum 
spin Hall system^- 11, 23[. Then, several new features of 
spin transport stemming from these helical edge modes 
have been also predicted lol- 13|. 

In Ref. 16j, we have clarified the ABS and tunneling 



near half filling [14J. The magnitude of spin-triplet /- 



wave pairing in this d x 2_ y 2 + /-wave pairing is enhanced 
by Rashba-type spin-orbit coupling originating from the 



conductance <jq in normal metal / NCS superconductor 
junctions for d xy +p-wave and d x 2_ y 2 + /-wave pairings. 
Both for them, new types of ABS appear, in stark con- 
trast to s + p-wave case. In particular, for d xy + p-wave 
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case, due to the existence of the Fermi surface splitting by 
spin-orbit coupling, there appears a Majorana edge state 
with flat dispersion preserving the time reversal symme- 
try. Reflecting the Majorana edge state, <jq has a zero 
bias conductance peak (ZBCP) in the presence of the 
spin orbit coupling. 

In the present paper, we study the time-reversal in- 
variant Majorana edge state in details. In particular, we 
examine the local density of state at surface, i.e. surface 
density of state (SDOS) for NCS superconductors with 
s + p, d x 2_ y 2 + /, and d xy + p-wave pair potential, re- 
spectively, based on the lattice model Hamiltonian. For 
d xy +p-wave pairing, we confirm that the existence of the 
special ABS appears as a zero energy state, due to the 
Fermi surface splitting by the spin-orbit coupling. The 
present ABS is a single Majorana edge mode preserving 
the time reversal symmetry. 

We also study the topological nature of the Majorana 
edge mode with flat dispersion. The ABS found here 
is topologically stable against a small deformation of 
the Hamiltonian if the deformation preserves the time- 
reversal invariance and the translation invariance along 
the direction parallel to the edge. We introduce a topo- 
logical invariant number ensuring the existence of the 
zero energy ABS, and clarify the relevance to the number 
of Majorana edge modes. It is revealed that the absolute 
value of the topological number equals to the number of 
the Majorana edge modes. 

Topological aspects of edge states have been attracting 
intensive interests in condensed matter physics. Espe- 
cially, it was highlighted by the discovery of the quantum 
Hall system (QHS) showing the accurate quantization of 
the Hall conductance (Jh which is related to the topolog- 
ical integer 0, [25| . It is known that chiral edge state 
is generated at the edge of the sample. The concept of 
the QHS has been generalized to the time-reversal (T) 
symmetric system, i.e., the quantum spin Hall system 
(QSHS) 



Greenfs function formalism. In Sec. Ill, the results of 
the numerical calculations of SDOS, topological invari- 
ant number, and SDOS in the presence of Zeeman mag- 
netic field are discussed. In Sec. IV, the conclusions and 
outlook are presented. 



FORMULATION 

In this paper, we consider the two dimensional square 
lattice with Rashba-type spin-orbit coupling. The model 
Hamiltonian is given by 
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where Ck a iyka) IS an annihilation (creation) operator 
for an electron, <r and <x; the Pauli matrices, and 
the energy dispersion of the electron on the square lat- 
tice, £k = — 2t(cosk x + cosfc y ) — (i with the nearest 
neighbor hopping t and the chemical potential /j,. The 
second term is the Rashba spin-orbit coupling, g(k) — 
(sink v , — sin k x , 0), and the third one is the pair poten- 
tial, Afc = iip(k)a y + id(k) ■ &<j y . In the presence of 
the Rashba spin-orbit coupling, the Fermi surfaces are 
split into two, and we suppose intraband pairings in each 
spin-split bands. Then the (i-vector of pairing function 
for triplet pairings d(k) is aligned with the polarization 
vector of the Rashba spin-orbit coupling, d(k) \\ g(k). As 
a result, the triplet component of the energy gap func- 
tion is given by d(k) = A t f(k)g(k) while that of singlet 
component reads ip{k) — A s /(fe). Here /(fe) is given 
by /(fe) = 1, sin k x sin k y and {cosk x — cosk y ) for s+p, 



2J. In QSHS, there exist the helical edge d xy +p and d^-^+Z-wave, respectively. [H We also in- 



modes, i.e., the time-reversal pair of right- and left-going 
one-dimensional modes. The above edge modes are gen- 
erated from non-trivial nature of bulk Hamiltonian and 
topologically protected. Furthermore, recently, to pursue 
analogous nontrivial edge state including Majorana edge 
mode in superconducting system has become a hot issue 
27]. Although there have been many studies about edge 
modes _m topologically non-trivial superconducting sys- 
the relation between edge modes (ABSs) 



troduce the Zeeman splitting term in an applied magnetic 
field fisH for later use. From the above Hamiltonian, we 
have the following retarded Green's function in the infi- 
nite system, 



{G (k.uj)}^ = (uj + i V )h x4 -H(k), 



with 



2cS 



tems 

and the surface density of states has not been fully clar- 
ified. 

In the presence of the single Majorana edge mode, the 
SDOS has an anomalous orientational dependence on the 
Zeeman magnetic field. We reveal that the SDOS with 
zero energy peak is robust against magnetic field in a 
certain applied direction. 

The organization of the present paper is as follows. In 
Sec. II, we introduce the Hamiltonian and the lattice 



H{k) = 



£,k Afc 
Al -tk 



(2) 



(3) 



where = £fc^2x2 + Ag(fc) • & — fi^H ■ <r with the 4x4 
(2 x 2) unit matrix 14x4(^2x2)- 

To calculate SDOS, we construct the Green's function 
in semi-infinite system. In the actual numerical calcula- 
tion, we use the periodic boundary condition along the 
x-direction with a sufficiently large size of mesh. To pre- 
pare the (100) surface at x = ao, we introduce vacuum 
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FIG. 1: (100) surface of square lattice. White and black cir- 
cles show the sites without and with potential V , respectively. 



layers on the right side of the (100) surface as shown in 
Fig. [T] Here, it is sufficient to introduce the two vacuum 
layers since there is no long range hopping or long range 
pairing over three lattice constant in the present model. 
On the vacuum layers, we add the following term to the 
Hamiltonian, 



(4) 



Xi—a\,a2 <y 



where rii a — c? ia Ci a is a number operator at site i with 
spin a, and V is the on-site potential. In the limit V — > 
oo, no electron exists on the vacuum layers. 

First, we switch on the potential V only at the site 
x = a\. The Green function G 1R (xi,x 2 ;k y ,uj) in this 
situation satisfies the following equation, 

G 1R (xi,x 2 ; k y ,uj) = G 0R (xi,x 2 ;k y ,u}) 

+G 0R {x 1 ,a 1 ;k y ,u)VT 3 G 1R (a 1 ,x 2 ;k y ,u>), (5) 

where G 0R (xi, x 2 ; k y , w) is the Fourier component 
G QR (k,uj) with respect to k x , 



G QR {x u x 2] k y ^) = -L^G°«(fc, w )e 



ik x {x\—X2 ) 



(6) 



Here N x is the number of ic-meshes, and T3 is the Pauli 
matrix in the particle-hole space. 



T3 



hx2 






-I 2 x2 



(7) 



Eq. ©, 

G 1R (xi,x 2 ;k y ,uj) = G 0R (xi,x 2 ;k y ,uj) 

-& R (x u av,k v ,u){& iR (ai,a 1 -,ky,u)}- 1 
xG 0S (ai,i2;t^). (8) 

Then, we switch on the potential V at the site x = a 2 
as well. For the Green's function G 2R (xi, x 2 ; k y , oS) in 
this case, we have the following equation, 



G {x u x 2 ;k yi uj) = G (x 1 ,x 2 ;k y ,uj) 

+G 1R (x u a 2 ;k y ,uj)Vf 3 G 1R (a 2 ,x 2 ;ky,uj), (9) 

Therefore, taking the V — >■ 00 limit, we obtain the 
Green's function in the semi-infinite system, 

G 2R (xi,x 2 ]k y ,uj) = G 1R {xi,x 2 \k y ,Lu) 

-G 1R (xi,a 2 ; k y ,uj){G 1R (a 2 ,a 2 ; fc a ,w)} _1 
xG 1R ia 2 ,x 2 ;k y ,uj). (10) 

From © and (JTUJ), we can calculate the SDOS p s (u), 
which is given by the local density of states at x = do, 



1 



^{G 2 a R (a ,a Q ;k y ,uj)}, (11) 



k v a=l,2 

while the DOS in the bulk is given by 
1 



Pb(w) 



N X N V ^ 5 

y k^kSy a=l,2 



Im{G^(fc, W )}, (12) 



In the V — > 00 limit, we have the following solution of 



where N x and N y are the number of meshes along the 
x and the y direction, respectively. In the calculation 
presented in the following, we choose N x — N y = 2 13 . 
We set t = Aq = 1 for the unit of energy. The number of 
electrons per unit cell is 0.3 (// ~ —2.40). To guarantee 
the convergence, we use rj = 0.03i for the infinitesimal 
imaginary parts in Eq. @. 



RESULTS 

In this section, we first focus on the local density of 
state at the surface, i.e., SDOS p s (w) for various par- 
ing symmetry. Since spin-singlet and spin-triplet com- 
ponents of pair potential mix in general, we introduce 
a parameter r s , which denotes the ratio of singlet com- 
ponent. The parameter r s (0 > r s > 1) is defined as 
A s = r s A and A t = (1 - r s )A . 

The SDOS for s + p-w&ve case is plotted in Fig. [5] 
In this case, the bulk DOS Pb(w) always has a [/-shaped 
gap structure as shown in solid lines. The magnitude of 
the gap is given by | A s - A t |. For A s > A t , the SDOS 
p s (oj) also has a [/-shaped gap structure as shown in the 
dashes lines in Fig. [U(a) and (b). On the other hand, for 
At > A s with r s = 0.2 and r s = 0, the resulting p s (u>) 
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has a residual value at ui = as shown in FigsJSJc) and 
(d), respectively. To show up the SDOS more clearly, 
we also show the angle resolved surface density of states 
(ARSDOS) -E a =i,2 l ™{G 2 a«{ao,a ;k y ,u)}. As shown 
in Fig. Ela), ARSDOS shows a full gap structure without 
any inner gap state for singlet dominant case (A s > At). 
On the other hand, as shown in Fig. [3jb) , the ARSDOS 
for triplet dominant case (A* > A s ) has two branches of 
Andreev bound state, which are dubbed as helical edge 
modes 111, |l2( . The presence of helical edge modes inside 
the bulk energy gap induces the residual SDOS inside the 
bulk energy gap. 



structures of the bulk energy gap. The corresponding 
p s (oj) (dashed line) also has a similar V^-shaped gap struc- 
ture for singlet dominant case (A s > A t ) with r s = 1 
and r s = 0.8. As shown in Fig|5]Ja), there is no inner gap 
state in ARSDOS, while the bulk energy gap has a strong 
k y - dependence. On the other hand, for triplet dominant 
case (At > A s ) with r s = 0.2 and r s — 0, p B (ui) has two 
additional peaks as shown in dashed lines of FigsHJc) and 
(d) as compared to the bulk density of states [16| (solid 
lines). The additional two peaks in the SDOS originates 
from the helical edge modes generated inside energy gap 
as shown in Fig. [5jb) (l6| . 
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FIG. 2: (Color online) Local density of states for s+p-wave 
in the bulk (solid lines) and at the surface (dashed lines) for 
A = 0.5 and r s = 1.0, 0.8, 0.2, and 0.0. 



FIG. 4: (Color online) Local density of states for d x 2_ y 2+f- 
wave in the bulk (solid lines) and at the surface (dashed lines) 
for A = 0.5 and r s = 1.0, 0.8, 0.2, and 0.0. 
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FIG. 3: (Color online) Angle resolved local density of state of 
s+p-wave pairing is plotted as a function of k y with A = 0.5. 
(a) A s = A and At = 0, (b)A t = A and A s = 0. 
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FIG. 5: (Color online) Angle resolved local density of state 
of d x 2_ y 2+ /-wave pairing is plotted as a function of k y with 
A = 0.5. (a)A s = A and A t = 0, (b)A t = A and A s = 0. 



Next we look at SDOS for d x 2_ y 2 + /-wave case plot- 
ted in Fig. SJ In this case, Pb(w) always has a V^-shaped 
gap structure as shown in solid lines reflecting the nodal 



The SDOS for d xy + p-wave NCS superconductors is 
plotted in Figs. |B] and [5] In this case, p s (u>) has very 
different line shapes as compared to former two cases. In 
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particular, for the spin-triplet dominant pairing, we find 
that the SDOS is very sensitive to the spin-orbit cou- 
pling. To show this, we first start with the case without 
spin-orbit coupling i.e. A = 0. See Fig. [5] pt,{u>) has a 
F-shaped gap structure as shown in solid lines, reflecting 
on the nodal structure of pair potential. The correspond- 
ing p s (uj) (dashed line) has a zero energy peak (ZEP) for 
singlet dominant case with r s = 1 and r s — 0.8. This 
ZEP comes from the mid gap Andreev bound state with 
flat dispersion as shown in Fig. [7(a) , and it is essentially 
the same as that appears in the surface state of high-T c 
cuprate O n the other hand, for triplet dominant 

case with r s — 0.2 and r s = 0, the ZEP disappear, but 
Ps{oj) supports two additional peaks instead, as shown in 
dashed lines of Figs[6(c) and (d), respectively. The ad- 
ditional two peaks are generated by the anomalous ABS. 
In the absence of Rashba spin-orbit coupling, the dis- 
persion of the anomalous ABS for d xy +p-wave pairing is 
very similar to the helical edge modes for d x 2_ y 2 +/-wave 
pairing, as shown in Fig. [7f b) . However, the intensity 
of ARSDOS near k y ~ is very low since the magni- 
tude of the excitation energy of the ABS is close to the 
bulk energy gap. In the 2D free electron model for NCS 
superconductor, it is shown that the dispersion of ABS 
corresponds the bulk energy gap and the intensity of AR- 
SDOS is completely absent for \k y \ < fc c . 16]. Thus, in 
contrast to the helical edge mode in the d x 2_ y 2 + /-wave 
pairing case, the values of the SDOS d xy +p-wave pairing 
at u/A = is very close to zero. 
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FIG. 7: (Color online) Angle resolved local density of state of 
d xy +p-wa,ve pairing is plotted as a function of k y with A = 0.0. 
(a) A s = A and At = 0, (b)A* = A and A s = 0. 



Ph(w) shows a F-shaped gap structure similar to Fig [5] 
Then, for singlet dominant case with r s = 1 and r s — 
0.8, the SDOS p s (uj) (dashed line in Figs. 15(a) and (b), 
respectively) has a ZEP similar to Figs. [6(a) and (b). On 
the other hand, for triplet dominant case with r s = 0.2 
and r s = 0, in addition to two peaks similar to those in 
Figs. [6(c) and [6(d), ZEP appears as shown in Figs[8(c) 
and[5(d). It is remarkable that ZEP is generated by spin- 
orbit coupling A [lij . 
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FIG. 8: (Color online) Local density of states for d xy +p-w&ve 
in the bulk (solid lines) and at the surface (dashed lines) for 
A = 0.5 and r s = 1.0, 0.8, 0.2, and 0.0. 



FIG. 6: (Color online) Local density of states for d xy +p-w&ve 
in the bulk (solid lines) and at the surface (dashed lines) for 
A = 0.0 and r a = 1.0, 0.8, 0.2, and 0.0. 



Let us now consider the d xy +p-wave case with nonzero 
spin-orbit coupling A. As shown in Fig. [5J the bulk DOS 



To show up the new ZEP much more clearly, we plot 
ARSDOS of 

d X y + p-wave case in Fig. [9] with A — 0.5. 
Comparing with the ARSDOS in FigO for triplet dom- 
inant case, we find that an additional zero energy state 
(ZES) appears. In the presence of the spin-orbit cou- 
pling, the Fermi surface is split into the large one with 
the Fermi momenta hi and the small one with k\. The 
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ZES exists only for k y between the split Fermi surfaces, 
namely, for k y with k% > \k y \ > k±. It can be shown 
that the Bogoliubov quasiparticle creation operator 7J 

for the ZES satisfies 7* = T-fe-[l3| Therefore, the ZES is 
identified as a Majorana Fermion. 

For singlet dominant d xy + p-wave NCS superconduc- 
tors, we obtain the SDOS and ARSDOS depicted in 
Figs[8ja)-(b) and Fig[9ja), respectively. At first sight, 
they look very similar to those in Figs|7ta)-(b) and FigJU 
However, for k y with ki > \k y \ > kx, we again have a 
single branch of a zero energy state on the edge. As is 
shown later, because of the existence of the time-reversal 
invariant Majorana Fermion (TRIMF), the SDOS for the 
d X y + p-wave NCS superconductor shows a peculiar de- 
pendence on the Zeeman magnetic field. 

Unlike the Majorana Fermions studied before, the 
present Majorana Fermion is realized with the time re- 
versal invariance. The TRIMF has the following three 
characteristics, (a) It has a unique flat dispersion. To 
be consistent with the time-reversal invariance, the sin- 
gle branch of ZES should be symmetric under k y — > —k y . 
Therefore, by taking into account the particle-hole sym- 
metry as well, the fiat dispersion is required. On the 
other hand, the conventional time-reversal breaking Ma- 
jorana edge state has a linear dispersion, (b) The spin- 
orbit coupling is indispensable for the existence of the 
TRIMF. Without the spin-orbit coupling, the TRIMF 
vanishes, (c) The TRIMF is topologically stable under 
small deformation of the Hamiltonian. The topological 
stability is ensured by a topological number, which will 
be shown below. 
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FIG. 9: (Color online) Angle resolved local density of state of 
dxy+p-wave pairing is plotted as a function of k y with A = 0.5. 
(a)A a = A and A t = 0, (b)A t = A and A 3 = 0. 



To see the topological nature of the TRIMF, we start 
from the Bogoliubov-de Gennes (BdG) Hamiltonian in 
the bulk system defined in cq. without Zeeman mag- 
netic field. In the Nambu representation, the BdG Hamil- 
tonian ([3]) has the particle-hole symmetry, 



CU{k)& = -U*{-k), C 



I 2x , 

hx2 



(13) 



In addition, from the time-reversal invariance, the BdG 
Hamiltonian satisfies 



(-W(fc)6 + = n*(-k), e= ( lcr Q v f ° 



(14) 



Therefore, we can define the operator T which anticom- 
mutes with the BdG Hamiltonian, 



{W(fc),f}+ = 0. 



(15) 



Here T is defined as the product of particle-hole trans- 
formation operator C and time reversal operator O, 



r = - -/("'(-) : : j °y 

<7y 



Now we take the basis which diagonalizes f 



U —^2x2 



with the unitary matrix Uy 



(16) 



(17) 



(18) 



Then, we find that the BdG Hamiltonian H(k) becomes 
off-diagonal in this basis, 



U r H(k)U T = 



q(k) 
g(fc)t 



(19) 



where q(k) = i k a y - A k . 

To classify the existence or nonexistence of the ZES 
with k y = ky, we change k x from — 7r to 7r for a fixed 
value of k y — ky. The winding number W is defined as 
a number of revolutions of det|g(fe)| = mi(fe) + im2(k) 
around the origin of complex plane when k x changes from 

— 7T to 7T, 



dO(k) 



dk x 



dk x , 



(20) 



with 6{k) = argdet|g(fc)| = tan- 1 (m 2 (fc)/mi(fc)). The 
resulting W(ky) must be an integer since the starting 
point and the end point of integration route are equiva- 
lent in the Brillouin zone. Therefore, the value of W(ky) 
changes discretely with the change of fc y -value. When the 
gap of the system closes at certain points on the integra- 
tion route, W(ky) is ill-defined since det|q(fe)| becomes to 
zero on such points. Thus, the value of W(k y ) does not 
change by the small change of fc y -value unless the energy 
gap closes. 

Here, we study the winding number for the d xy + p- 
wave NCS superconductors for on the two trajectories 
for ky — 0.2U7T and 0.457T as shown in Fig. [TU] For 
ky = 0.27T, the trajectory crosses the Fermi surfaces four 
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times. On the other hand, for k y = 0.457T, it crosses 
two times. As was explained above, the winding number 
W(k y ) does not change unless the energy gap closes. For 
d xy + p-wave parings, the energy gap closes at (±fci,0), 
(±fc 2 ,0), (0,±fci) and (0,±fc 2 ). Therefore, W(k y ) can 
change only at k y — ±fci, k y = ±fc2 or k y = 0. Thus, it 
is sufficient to study the representative point in each k y 
range. 



it 



A,0 




-71 71 



FIG. 10: Fermi surfaces of spin-split bands due to the Rashba 
spin-orbit coupling (solid lines) and cutting lines k y = 0.207T 
and 0.457T (dashed lines). 




k x 

FIG. 11: Polar plot of 9(k) for A s = A and At = as 
a function of k x with fixed k y , (a) k y — 0.2n and (c) k y = 
0.457T, where white circle corresponds to k x = ±7r and arrows 
show the direction in which the values of k x increase. The 
corresponding 9 — tan -1 (7712/7711) is plotted as a function of 
k x for (b) k y = 0.2n and (d) k y = 0.457T. 



We first focus on the singlet dominant d xy + p-wave 
NCS superconductor with A s — Aq and A t — 0. For 
k y = 0.27T, as shown in Fig. Illf a). mi and 7712 draw 
a curve which turns anti-clockwise twice around the 
(mi, vri2) = (0, 0) if we change k x from — tt to tt. At the 
same time, 9(k) changes twice from — tt to tt as shown in 
Fig. Illf b). Therefore, the resulting winding number W 
is W — 2. On the other hand, for k y = 0.457T, mi and m-i 
draw a curve which turns anti-clockwise once around the 
(mi, m 2 ) = (0, 0) as shown in Fig. [TTTc). Moreover, 9(k) 
changes once from — tt to tt as shown in Fig. Illf d). Thus 
the resulting winding number W equals to W = 1. This 
means that these two cases belong to different topological 
classes. In a similar manner, it is possible to generalize 
above argument for other k y a (—tt < k y < tt). We sum- 
marize the obtained results in Table fl] (a) . 

The same plot of rrii(k) and 9{k) for the triplet dom- 
inant d xy + p-wave NCS superconductors with A t = Ao 
and A s = is shown in Fig. [T^l For k y = 0.2tt, mi and 
mi draw a curve as shown in Fig. I12f a). However, it does 
not turn around the (mi, 777,2) = (0,0) with the change 
of k x from — tt to tt. At the same time, 8(k) does not 
change from — tt to tt as shown in Fig. I12f b). in contrast 
to that in Fig. ITTT b). This corresponds to the fact that 
the resulting winding number W equals to W = 0. On 
the other hand, for k y — 0.457T, mi and 7772 draw a curve 
which turns clockwise once around the (mi, 7772) = (0, 0) 
as shown in Fig. IT2T c). Also, 9{k) changes once from tt 
to — tt as shown in Fig. [T2T d). Therefore, the correspond- 
ing winding number W equals to W = — 1. In a manner 



to the singlet dominant case, it is possible to generalize 
this argument for other k y s. (—tt < k y < 7r)The obtained 
results are summarized in Table H^b). 

Comparing the obtained winding number W(k y ) in Ta- 
ble U with the ARSDOS in Fig [HI we notice that zero en- 
ergy ABS appears only for k y with nonzero W(k y ). This 
correspondence implies that the existence and the sta- 
bility of the zero energy ABS is ensured by the winding 
number. Indeed, we can say the absolute value of W , i.e., 
W I equals to the number of branch of zero energy ABS 



3l|. In other words, the number of TRIMF equals to the 
absolute value of W . Furthermore, the TRIMF found 
here are topologically stable against a small deformation 
of the BdG Hamiltonian if the deformation preserves the 
time-reversal invariance and the translation invariancc 
along the direction parallel to the edge, both of which 
are necessary to define the winding number. It should 
be remarked that both for the singlet dominant case and 
the triplet dominant one, a single TRIMS is generated 
for fc 2 >| k y |> ki. 

As is shown above, the time-reversal invariance is es- 
sential to define the winding number W ensuring the 
topological stability of the TRIMF. Therefore, in general, 
if we apply a perturbation breaking the time-reversal in- 
variance, the winding number W becomes the meaning- 
less. In other words, a time-reversal breaking perturba- 
tion may change the SDOS of the TRIMF substantially. 

As a time-reversal breaking perturbation, we consider 
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(a) 



ky 


k y > ki 


&2 > ky > fel 


fel > ky > 


> ky > -fel 


— fei > k y > —k-2 


— k2 > k y 


W(ky) 





1 


2 


-2 


-1 





W{ky) 





-1 








1 






TABLE I: The winding number W(k y ) for d xy + p-wave with (a) A s > A t , (b) A t < A s 



(a) k =0.20x 



(b) fc=0.20jr 
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m^k) 



FIG. 12: Polar plot of 0(fc) for A t = A and A s = as 
a function of k x with fixed fc H , (a) k y — 0.2n and (c) k y — 
0.457T, where white circle corresponds to k x = ±7r and arrows 
show the direction in which the values of k x increase. The 
corresponding 6 — tan _1 (m2/mi) is plotted as a function of 
k x for (b) k y = 0.27T and (d) k y = 0.457T. 



the Zeeman magnetic field. Let us look at the SDOS 
in the presence of Zeeman magnetic field H. First, we 
consider the singlet dominant d xy + p-wave NCS super- 
conductor. Without the spin orbit coupling, the Andreev 
bound state becomes conventional and can be expressed 
by double TRIMFs. In this case, as we expected, it is 
found that the zero energy peak of SDOS from the dou- 
ble TRIMF is split into two by any Zeeman magnetic 
field H as shown in Fig. [13] [32]. It is also noted that 
the resulting SDOS is independent of the direction of H 
due to the spin-rotational symmetry in the system. Thus, 
the SDOSs Ps (uj) for p B H x = O.lt (H || x), p B H y = O.lt 
(H || y) and HbH z — O.lt {H \\ z) are the same. 

On the other hand, in the presence of the spin-orbit 
coupling, this property is not satisfied anymore. As 
shown in FigfTH p s (w) has a different values for pbH x = 
O.lt [Fig. p B H y = O.lt [Fig. [He)] and p B H z = 

O.lt [Fig. Hlf d)] . respectively. The orientational depen- 
dence of SDOS is due to the presence of the spin-orbit 
coupling since the spin-rotational symmetry is broken. It 
is noted, three peak structure including zero energy peak 
appears for {1bH v = O.lt (H \\ y). The presence of the 



1 

0.5 



0, 



(a)/i B H=0 



(b)M B H=0.1t 




0.4-0.2 0.2 0.4-0.4-0.2 0.2 0.4 

co/A co/A 

FIG. 13: Local density of states at the surface for d xy +p- 
wave with A s = Ao and At = for A = and (a) /j,bH — 0, 
(b) ft B H w = 0.1 (H || x). 



remaining ZEP under the Zeeman magnetic field along 
the edge, which is H y in this case, is relevant to the ex- 
istence of the single Majorana fermion at &2 > \k y \ > k±. 
On the other hand, the double TRIMFs at \k y \ < ki is 
split into two for any direction of Zeeman magnetic field. 
Therefore, the height of ZEP p s (0) drastically decrease 
even if ZEP remains for pbH v = O.lt (H \\ y). 



(a)M B H=0 




(b)M B H x =0.1t 



1 

0.5 



0, 



(c)// B //=0.1f 



(d) ju B H z =0.1t 



-0.4-0.2 0.2 0.4-0.4-0.2 0.2 0.4 



CO/ A n 



CO/ An 



FIG. 14: Local density of states at the surface for d xy + p- 
wave with A s = Ao and At = for A = 0.5 and (a) pbH = 0, 
(b) p B H x = O.lt (H || as), (c) p B H y = O.lt (H \\ y), (d) 
psH z =0.1* (H || z). 



Next we consider the triplet dominant d xy + p-wave 
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NCS superconductor. For simplicity, we suppose that 
only the triplet component A p exists. As was shown in 
FigEJb), in the absence of the Zeeman magnetic field, a 
zero energy ABS exists and it can be described by a single 
TRIMF. The resulting p s (uj) has a sharp ZEP without H 
[Fig. I15f a)]. In a manner similar to the singlet dominant 
case above, the SDOS p s {u) under the Zeeman magnetic 
field has a strong orientational dependence of H as shown 
in Figs. [15] (b), (c) and (d). It is noted that ZEP remains 
when the applied Zeeman field is along y-direction. This 
implies that the single TRIMF in the d xy + p-wave NCS 
superconductor is robust against the Zeeman magnetic 
field applied in the direction along the edge. 



(a)// B ff=0 




(b)M B H x =0.lt 




(d) jU B H z =Q.lt 




-0.4-0.2 0.2 0.4-0.4-0.2 0.2 0.4 

co/A co/A 

FIG. 15: Local density of states at the surface for d^+p-wave 
with A t = Ao and A s = for A = 0.5 and (a) psH = 0, 
(b) fi B H x = O.lt (H || as), (c) fj, B Hy = O.lt (if || y), (d) 
PbHz = O.lt (H || z). 



We would like to emphasis that, as shown from these 
calculation, the tunneling spectroscopy with the Zeeman 
magnetic field is available to identify the single TRIMF 

m til 6 dxy 

+ p-wave NCS superconductor. Simultane- 
ous existence of the strong orientational dependence of 
magnetic field and the robust zero energy peak under a 
certain direction of the magnetic field is a strong evidence 
to support single Majorana Fermion in this material. 



SUMMARY 

In the present paper, we have obtained SDOS, edge 
modes of NCS superconductors by choosing s + p, 
d x 2_ y 2 + /, and d xy + p-wave pair potential based on 
the lattice model Hamiltonian. For d xy + p-wave pair- 
ing, special ABS appears as a zero energy state for 



&2 >| k y |> ki with wavenumber parallel to the inter- 
face fcy, where k\ and ki denote the magnitude of the 
Fermi wavenumber in the presence of spin-orbit coupling. 
The present ABS is a single Majorana edge mode pre- 
serving the time reversal symmetry. We have defined a 
new type of topological invariant number and clarify the 
relevance to the number of the time-reversal invariant 
Majorana edge modes. We have found that the abso- 
lute value of the topological number equals to the num- 
ber of the Majorana edge modes. The single Majorana 
edge mode is generally induced by the spin orbit cou- 
pling. In the presence of the single Majorana edge mode, 
the SDOS has a strong orientational dependence of the 
magnetic field. In a certain direction of the applied mag- 
netic field, the single Majorana edge mode is robust and 
the resulting zero energy peak of SDOS remains. The 
present future may serve as a guide to detect Majorana 
Fermion in non-centrosymmetric superconducting super- 
conductors by tunneling spectroscopy. 

This work is partly supported by the Sumitomo Foun- 
dation (M.S.) and the Grant-in-Aids for Scientific Re- 
search No. 22103005 (Y.T. and M.S.), No. 20654030 
(Y.T.) and No.22540383 (M.S.). 
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